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This paper presents a Stand-alone Hybrid Renewable Energy System 
(SHRES) as an alternative to fossil fuel based generators. The Photovoltaic 
(PV) panels and wind turbines (WT) are designed for the Malaysian low 
wind speed conditions with battery Energy Storage (BES) to provide electric 
power to the load. The appropriate sizing of each component was 
accomplished using Non-dominated Sorting Genetic Algorithm (NSGA-II) 
and Multi-Objective Particle Swarm Optimization (MOPSO) techniques. The 
optimized hybrid system was examined in MATLAB using two case studies 
to find the optimum number of PV panels, wind turbines system and BES 
that minimizes the Loss of Power Supply Probability (LPSP) and Cost of 
Energy (COE). The hybrid power system was connected to the AC bus to 
investigate the system performance in supplying a rural settlement. Real 
weather data at the location of interest was utilized in this paper. The results 
obtained from the two scenarios were used to compare the suitability of the 
NSGA-II and MOPSO methods. The NSGA-II method is shown to be more 
accurate whereas the MOPSO method is faster in executing the optimization. 
Hence, both these methods can be used for techno-economic optimization of 
SHRES. 
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1. INTRODUCTION 

In recent years, the increasing concern on the depletion of fossil fuel and global warming has 

catalyzed the growth of renewable energy sources due to their promising economic and environmental 
benefits [1], [2]. Wind turbines and solar photovoltaic are commonly used in the renewable energy system to 

supply power to consumers in the remote regions because there is no fuel cost involved, easy to install and 

are also non-polluting. Nevertheless, designing a renewable energy system can be a challenge. Thus, 
knowledge of all aspects that influences system performance and component sizing is a precondition for an 
accurate SHRES design. Large fluctuations in climatic and meteorological conditions cause intermittency of 
renewable energy sources. Malaysia which lies close to the equator has seasonal wind speed and does not 
have a comprehensive wind assessment [3]. This problem can practically be overcome by using battery 
energy storage system [4]. They design and optimization methods are important aspects for SHRES to 


Journal homepage: http://iaescore.com/journals/index.php/IJPEDS 








464 □ 


ISSN: 2088-8694 


guarantee supply reliability and security, and also to ensure maximum utilization of PV panels, wind turbines 
and battery energy storage, based on the load profile [5], [6]. 

There are several methodologies that utilize traditional optimization methods to design a techno- 
economic hybrid energy system based on the Loss of Power Supply Probability (LPSP) [7]-[9]. A common 
drawback of these optimization methods is the low calculation efficiency, therefore consuming excessive 
central processing unit time. In addition, the optimization methods cannot find the best compromise point 
between the objective functions. On the other hand, Artificial Intelligence (AI) methods are able to achieve 
all the conditions, such as LPSP and COE. Authors in [10]-[12] presented a methodology to calculate the 
optimal number of PV panels, wind generation, and battery using genetic algorithm approach by calculating 
LPSP and system cost. In [13], AI methods such as the Non-dominated Sorting Genetic Algorithm (NSGA- 
II) and Multi-Objective Particle Swarm Optimization (MOPSO) were used to produce a Pareto-optimal 
solution in a single simulation run. Several authors have studied the hybrid optimization system using 
NSGA-II. In 2017, Moslem Yousefi et al. [14] used NSGA-II algorithm and HOMER software to find the 
robust project design of SHRES by the varying engine loads with optimal Annual Energy Recovery (AER) 
and total cost of the system. Much attention was paid to discuss SHRES and Battery Energy Storage (BES) 
sizing. Reference [15] discussed the economic approach of multi-optimization of a standalone hybrid PV- 
Wind - Battery and diesel generator system through the application of multi-objective using NSGA II 
method. The economic benefits include the minimization of power generation cost and maximizing the useful 
life of the battery, including the life loss, fuel, environmental, and maintenance cost. In addition, it has 
considered the lifetime characteristics of lead-acid batteries. 

Ce Shang et al. [16] focused on the battery energy storage system sizing in stand-alone hybrid 
power system to guarantee reliability and minimize the levelized cost of energy using NSGA II method. In 
reference [17], techno-economical optimization for HRES was applied using NSGAII method to analyze the 
trade-off between three conflicting objectives: total cost, autonomy level, and wasted energy rate. However, 
the optimal sizing of the system components was not considered. The hybrid solar/wind system with the 
traditional fossil fuel-fired generators was described in [18]. While two objectives that controlled the NSGA¬ 
II procedure was proposed to minimize the cost and emission. To achieve the best compromise solution, the 
Fuzzy priority ranking has used. The paper presented effectuality of the algorithm for evaluating through 
solving cost and emission dispatch issue without considering the power generation, for comparison reasons 
and results were compared with methods contained in the literature. Reference [19] presented small hybrid 
renewable system depends on the cost and environmental criteria using NSGA II technique, with two 
integrated energy storage units from battery banks and hydrogen storage system combined. Hence, 
minimized the COE and greenhouse gas emission (C02).The main contribution of this work was that the 
computing of total greenhouse gas emissions according to life cycle analysis of each system’s component. 

Many researchers have investigated the hybrid system using MOPSO method. Authors in [20] 
presented the optimization of an off-grid hybrid micro-grid system to determine the optimal sizing in twelve 
Swedish regions. The optimal design was selected after running the multi-objective optimization method to 
determine the trade-off between the three objectives: LPSP, COE, and the environmental impact (C02). A 
hybrid renewable energy system with multi-storage system configuration for buildings in Canada that adopts 
the MOPSO method was proposed in [21] for optimal economic operation in order to minimize the total Net 
Present Cost (NPC) and C02. Reference [22] proposed the design of a stand-alone hybrid generating system 
to determine the optimum sizing of the number and type of PV panels, wind turbines, battery bank, as well as 
diesel generators location using MOPSO method. The sizing was done based on a one-year data to minimize 
the cost and emission. 

The NSGA-II and MOPSO are the modem random optimization methods that are able to find 
Pareto. Hence, both these methods are applied to design the SHRES and to minimize the multi-objectives 
such as LPSP and COE. This paper presents a comparative evaluation of the performance of NSGA-II and 
MOPSO to determine the optimal sizing of SHRES using Pareto optimization. In order to determine the best 
combination of energy sources and to ensure their seamless integration into the distribution system to be at 
the optimal size, the numbers of PV panels, WT system, and batteries are used as the decision variables to 
minimize the LPSP and COE. The rest of the paper is organized as follows. A brief overview of hybrid 
WT/PV models and BES is presented in Section 2. The optimal configuration of the stand-alone hybrid 
system is explained in Section 3. The comparative analyses of the simulation are discussed in Section 4, 
followed by the conclusion in Section 5. 


2. POWER CIRCUIT OF THE (SHRES) 

In order to predict the SHRES performance, the energy sources need to be practically designed to 
meet the load demand. At the same time, the power obtainable from a hybrid renewable system has 
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significant fluctuations due to weather conditions and hence, the constant load demand may not be met. To 
mitigate this issue, a battery bank can be integrated to the hybrid system. However, the high cost of batteries 
is an issue in renewable energy systems. Thus, optimizing the size of the PV-WT-BES system becomes 
essential. These contributions reduce the capital cost and increase the chances of investment in renewable 
energy plant installation. Therefore, the optimal combination of renewable power resources with appropriate 
storage sizing, as proposed in this work, will give a vital contribution for the future economic feasibility of 
such plants, thus making the design more attractive for investors. 

2.1. Wind turbine modeling 

The wind is characterized by its speed and direction and is affected by factors, such as geographic 
position, meteorological factors and height above ground level. Wind turbine reacts to the wind, capturing a 
part of its kinetic energy and switching it into usable energy. The output power of wind turbine is determined 
as a function of the rated wind speed (V T ), the cut-in wind speed (V C i) and the cut-out wind speed (V co ) 
according to the following (1): 

f 0 V < V ci or V > V co 

p "- ^ivsv r at 

l P r Vr < V < V co 

Where, Pwt is the output power by wind turbine, P r is the rated wind power, V is the wind 
speed, V ci> V r and V co represent the cut-in wind speed, nominal wind speed, and cut-out wind speed 
respectively. The turbine cut-in speed is small, which enhances the effective operation of the system even 
under low wind speed [20]-[23]. 

2.2. Solar PV array modeling 

Solar panels are defined as a group of cells connected in parallel and series to generate the required 
electrical power based on meteorological factors such as solar radiation and temperature. The current model 
used to predict the output power of a PV module can be expressed through the following (2) and (3) [24]: 


Ppv= PVstc (ffj (l + K t (T c - T ref ) ). 

(2) 

T c = T a + (0.0256 * G). 

(3) 


Where, PVstc is the nominal power in (kW), G is the global solar radiation (kW/m2), G re f is the 
solar radiation under STC (1000/m 2 ), Tc is the temperature of PV cell, T re f = 25 C°, K T is the PV temperature 
coefficient, 3.7*10 -3 (1/C) and T a is the surrounding temperature. 

2.3. Battery storage modeling 

The typical batteries that are used for hybrid energy system in areas of low wind speed and 
intermittent solar radiation conditions are the lead-acid and lithium-based batteries[25]. Commonly, both 
these batteries are employed in most large-scale energy storage projects because of their low cost, long life 
span, and durability, in addition to their commercial availability [26]. Membrane based lead acid batteries are 
also available in the market presently. Battery storage is sized to meet the load demand during a shortage or 
interruption of renewable energy source, usually referred to as Autonomous Days (AD). Regularly, AD is 
taken to be one to three days according to (4). Typically battery capacity design depends on the load and AD. 
Thus, battery capacity can be calculated using the following equation: 


C B = (E l *AD) / (DOD* p bat * p inv ) 4) 

Where, Cb is the battery capacity, El is the load demand, DOD is the depth of discharge, gbat is the 
efficiency of battery and jii nv is the inverter efficiency. It is noteworthy that, the preceding expression is only 
used when the hybrid PV/WT system is unable to supply the required energy [20]. The hybrid PV/WT and 
battery system is formulated as a multi-objective optimization problem to improve the techno-economic 
performance simultaneously. 
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2.4. Hybrid energy management system 

The uncertainty of renewable energy supplies (Ere) is the main limitation of hybrid renewable 
energy plants. Therefore, an energy management strategy is required to complement the exchange of power 
from the generating sources and the load under variable weather conditions. It is calculated using the 
following (5) and (6): 


AEnet (t) = Ere (t) - El (t) (5) 

Ere (1) = Nwt* Ewt (t) + Npy * Epy (t) (6) 

Where, Apnet is the net energy of SHRES, Nwt is the number of wind turbines, Npv is the number of 
PV panel, E W t is the energy generated by the wind turbines, Epy is the energy generated from the PV panels, 
and E l (t) is the load demand at hour (t) where (t) equals one hour. 

The following cases are taken into account in this article, to simulate an energy management strategy, as 
depicted in Figure. 1: 

a. When the generated power is higher than the load demand, the surplus power is employed to charge the 
battery bank. 

b. When the generated power is higher than the load demand and the state of charge of the battery bank is 
full, the surplus energy is consumed in a dump load. 

c. When the generated energy is lower than the load demand, the battery bank is discharged to sufficiently 
supply the load demand. 



Figure 1. Flowchart of the hybrid energy system 


3. OPTIMAL CONFIGURATION OF THE STAND-ALONE HYBRID SYSTEM 

Once the hybrid component specifications have been determined, two cases will be investigated 
based on multi-objective optimization using NSGA-II and MOPSO methods. In the first case, the hybrid 
system will consist of PV panels, WT, and battery bank, while the second case consists of PV panels and 
battery bank only. The following two sub-sections illustrate the definitions of the objective functions 
in detail. 

3.1. Reliability analysis 

There are two approaches to determine the long-term performance of LPSP in a stand-alone hybrid 
system, namely, chronological method and probabilistic techniques [27]. The chronological method is more 
common and accurate, especially to determine the energy produced from the battery and the computational 
time is typically larger than that of probabilistic models. It is common in the chronological models to perform 
a one-year simulation with a one-hour time step. Computation time is especially necessary because this kind 
of model is generally used for component size optimization that requires several iterations. Hence, the 
chronological method is utilized in this paper. The LPSP is defined as the probability of unmet load over the 
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total energy produced [26], as mentioned in the first objective. The unmet load can be calculated by utilizing 
the deficit power between the load and sources in SHRES through the following (7) [3] : 

LPSP= ^ t=l E L (0 ~ e RE (0 + ( Eb (t) -Ebmin ( 

iLl *l(0 v ’ 


3.2. Economic analysis 

COE is defined as the average cost per kilowatt-hour ($/kWh) of electrical energy produced by the 
hybrid energy system [24], which can be achieved via the following equations: 


COE = (CRF * TAC) / E L 


( 8 ) 


CRF = 


d (l+i) T 
(l+i) T -l 


(9) 


Where, CRF is the capital recovery factor, which calculates the present value of system components 
by considering the interest rates (i = approximated as 6%) and project life span (T= 20 years in our case). 
TAC is the total annualized cost in $. The total annualized cost is the sum of the annualized capital cost (CC), 
operation and maintenance cost (Co&m) and replacement cost (CR). 


TAC = C c + Co&m + Cr (10) 

Cc = (Npv*Ppv*Cpv) + ((Nwt*Pwt*Cwt) + (Cwt*Nwt* 20/100)) + (Ns*Cb*CB) + 

(Ninv * Cinv) + (Cregpv + Creg_wt) • (11) 

Where, Cpy, Cwt, and Q>, are PV panel unit-price, WT unit- price, and battery unit-price 
respectively. Also, Cinv, Creg_pv, and Creg_wt are inverter unit price, regulator of PV price, and regulator of 
WT price respectively. In addition, Ninv is the number of inverters and equal 2 as well as the cost of the wind 
tower estimated as 20% of the system capital cost. The Co&m taken as 1% of the total cost as earlier supported 
[28]. For the Replacement Cost, expressed in (16) and (17). 


CR = CREP*SFF(i,P R _LF) (12) 

Crep = i r * ((Nb * Crb) + (Ninv * Cr.inv) + (Nregpv * Cr_reg_pv) 

+ (Nreg wt * Cr_reg_wt)) (13) 

Where SFF expressed as the sinking fund factor, Pr lf is lifespan of components (battery, inverter, 
and regulators). Crb, Crinv, Cr_reg_pv, and Cr_reg-wt are the replacement cost of the battery, inverter, PV 
regulator and WT regulator respectively. 

While Nreg pv is the number of voltage regulator and equal 1, Nreg_wt is the number of WT 
regulator and equal 1. Whereas the sinking fund factor is clarified by the following equation: 

Therefore, the sinking fund factor is clarified by the following equation: 


SFF (i, Pr lf) = (1+ijP R LF _ 1 (14) 

3.3. Solution methodology 

MOPSO and NSGA-II algorithms work by creating new reliable solutions. However, their working 
mechanism differs. Table! shows the differences between the two algorithms. 


Table 1. Difference between NSGA-II and MOPSO algorithms 
NSGA II MOPSO 


Selection, crossover, and mutation are used during each 
generation. Those individuals or chromosomes are combined to 
create children model. 

Great at finding the global optimum solution. 

More complex due to mutation and crossover; takes extra time 
compared to MOPSO. _ 


Particle positions are affected by their self-data and information 
sharing among swarm member. It uses two equations: velocity 
and position. 

Capable of finding the local optimum 

Fast and easy to implement as only a few parameters need 
adjustment. _ 
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In recent years, the implementation of NSGA-II and MOPSO algorithms based on GA and PSO are the 
foremost techniques adopted for global optimization in various fields, such as business, engineering, and in 
stochastic nature of renewable energy applications [1]. In this study, the NSGA-II and MOPSO methods are 
utilized to size the stand-alone hybrid PV/wind systems and these methods are shown to be better than the 
single-objective methods such as hybrid genetic algorithm (GA) or particle swarm optimization (PSO). 

3.3.1. Optimal configuration based on NSGA II algorithm 

In 2002, Deb proposed the use of NSGA-II algorithm[29], in which the population is distributed into 
several non-domination levels and each solution is assigned a fitness equal to its non-domination level. The 
algorithm can be summarized as follows [19]: 

a. In the first step, the requested input data are provided. This data involves the specifications of the hybrid 
system (hourly radiation, temperature and wind speed) and also load demand, data to compute the 
technical and economic functions and data to assess the constraint situations. 

b. Upper and lower bound of the number of PV-WT and BES 

c. The energy output of PV and wind turbine are calculated through the PV and wind models by using (1- 
3). The model of energy storage battery by using (4) with the total capacity (CB) is allowed to charge 
and discharge up to a limit defined by the maximum depth of discharge (DOD), by using (5, 6) and 
Figure 1. 

d. A random parent population (Pi) is created with size N. Then, the population of children (Qi) including 
N solutions is produced through genetic manipulation (crossover and mutation). 

e. Calculate the objective functions for each individual of Pi population (LPSP and COE) using (7-14). 

f. The two populations are combined to form the (Ri) population with size 2N. 

g. Classification of the Ri population is made in accordance with the Pareto front on the bases of fitness 

(non-dominated sorting is performed to determine the rank (front) of each population member). 

h. The next population of one of the fronts is built according to priorities by performing a general 
comparison of the members of the Ri population. 

i. Since the size of Ri is equal to 2N, the remaining solutions can simply be ignored because it is 
impossible to place all members in the new population (Pi+1). 

j. The end, if the maximum number of iterations referred in Step 3 is reached, the non-dominated sorting 

resolution at the last iteration was achieved as the optimal sizing and design for the SHRES. Otherwise 

go back to number 2. 

Figure 2 displays the NSGA-II algorithm application approach. This approach is aimed at finding 
the numbers of variables that depend on the state of Npv, Nwt, and Nbes in order to meet the minimized 
system LPSP and COE. This is executed by using NSGA-II optimization toolbox in MATLAB, with a 
selected population size of 200, the crossover value of 0.8, and the maximum simulation generation number 
set at 500. 



Figure 2. (a) and (b) NSGA-II procedure [27] 
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3.3.2. Optimal configuration based on MOPSO algorithm 

In 1995, Kenney and Eberhart showed that the Particle Swarm Optimization (PSO) has two separate 
concepts: a) social interaction which is exhibited by swarming and b) field of evolutionary calculation. In 
PSO, the two best values will determine the position of each particle. 

PSO is an AI approach which is based on the swarm social interaction within a field of evolutionary 
calculation, as proposed [30]-[32]. The algorithm determines the two best positions for each particle. 
Initially, the best value is obtained, called the individual best (pbest), and is retained by the particle, while the 
next value is determined by the PSO optimization algorithm within a global best populations. Individual 
particle position defines the particles variable target values and velocity that is applied to monitor the overall 
global best value (gbest). The fitness equation of this algorithm is to search the best solution from amongst all 
possible available options, with additional constraints added. The algorithm is based on each particle fitness 
appraisal, individual and global best fitness update, alongside with particle position and velocity. 

During the operation of the algorithm, each particle keeps the best fitness value that it has achieved. 
The particle with the best fitness value is calculated and updated during iterations. In this case, each particle 
represents a potential configuration of the PV-wind turbine and battery hybrid system: Npv, Nwt and Nbes, 
and the search space dimension are three. Then, the objective function of each particle is computed, 
corresponding to each scenario configuration (LPSP and COE). The following steps illustrate utilized this 
method for TIRES as the following: - 

a. In the first step: Initialization, the requested input data are provided. This data involves the 
specifications of the hybrid system (hourly radiation, temperature and wind speed) and also load 
demand, data to compute the technical and economic functions and data to assess the constraint 
situations. 

b. Upper and lower bound of the number of PV-WT and BES 

c. The energy output of PV and wind turbine are calculated through the PV and wind models by using 
(1-3). The model of energy storage battery by using (4) with the total capacity (CB) is allowed to charge 
and discharge up to a limit defined by the maximum depth of discharge (DOD), by using (5, 6) 
and Figure 1. 

d. Constants: 

Personal and global coefficients, Cl = C2 =2. 

Inertia weight, w = 0.9. 

e. The position and velocity of particles are randomly selected in order to generate the initial population 
and then applied to the objective functions to find the optimum fitness value, by using (7-14). 

f. Evaluate the fitness value, with minimum LPSP and COE 

g. Calculate and update (Pbest and gbest) 

h. Calculate and update velocity and position of each particle 

i. Apply the updated value to find optimum value of LPSP and COE 

j. The end, if the maximum number of iterations referred in Step 3 is reached, the non-dominated sorting 
resolution at the last iteration was achieved as the optimal sizing and design for the SHRES. Otherwise 
go back to number 2. 

The sizing model of the hybrid PV/wind energy systems is more complex than the single-source 
generation systems. This is because the variables must be considered for system optimization. Moreover, 
system performance over a long-term, economic parameters, and reliability objectives should be well 
thought-out in order to achieve a suitable compromise between COE and LPSP. NSGA-II and MOPSO are 
the appropriate methods with regards to global optimization and the random nature of renewable power 
sources. These methods have been used in many hybrid applications in recent years [4], [33], [34] 

In this paper, the NSGA-II and MOPSO algorithms were applied to size the stand-alone hybrid 
PV/WT and battery system. In principle, these algorithms aim to find the optimum number of PV panels, 
wind generation, and batteries to minimize the LPSP and COE. The optimization method was done using 
MATLAB software by executing the same number of iterations and population in both algorithms. Finally, 
the results were compared and the algorithm with the best result was identified. 

The optimization method starts with the following input values: hourly output power for PV, WT 
and load profile. Then, the loop iterates the NSGA-II or MOPSO algorithms. The optimum results generated 
by the two algorithms are converted to the nearest integer values. This gives the sizing of the generating 
source. The optimum sizing along with COE and LPSP values are recorded for each run and stored in an 
array for each iteration. Finally, all appropriate solutions are obtained through non-dominated optimal called 
Pareto front [35]. The number of each generating units are provided. Flow charts of the NSGA-II and 
MOPSO algorithms are shown in Figure. 3. 
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Figure 3. NSGAII or MOPSO optimization flow chart for stand-alone hybrid energy system 


4. RESULT AND DISCUSSION 

The methodology was applied to find the optimal size for SPIRES. The load profile for a typical 
rural village in Malaysia consisting of 20 households is shown in Figure 4 and the total energy consumption 
per day is 138.4 kWh [36]. The maximum solar radiation was approximately 1050 W/m 2 and the maximum 
wind speed was recorded at 5m/sec. The annual meteorological conditions in Malaysia, and solar radiation 
and wind speed are illustrated in Figure 5a and 5b, respectively. The data was obtained from the Malaysian 
Meteorological Department. Table 2 indicates the parameters of PV panels, wind, battery and inverters. 
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Figure 4. Hourly load profile for rural area in malaysia 
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Figure 5. Annual meteorological conditions for rural area in malaysia, 
(a) Solar radiation and (b) wind speed 


Table2. Parameters of PV, wind turbine, battery and inverter 

PV modules 

Specifications 

Wind turbine 
specifications 

Battery 

specifications 

Inverter specifications 

Power max = 320W 

Rated voltage = Vmpp 54.7v 
Rated current = Impp 5.49A 
Initial cost = $ 290 [37, 38] 
PV regulator cost =$ 750 [39] 
Life time = 20 years 

Rated output power = 3 kW 
Generator voltage = 230 V-ac 
Cut-in wind speed = 2 m/s 

Rated wind speed = 12 m/s 

Initial cost =$ 2800 [40] 

Wind regulator cost = $ 750 [41] 
Life time = 20 years 

Rated capacity = 

1000Ah 

Rated voltage =2V 
Efficiency = 85% 

DOD = 70% 

Initial cost =$ 230 [42] 
Life Time =10 Years 

Rated output power 6500 W 
Input Voltage 12Vdc/ 24Vdc 
Frequency 50 HZ 

Efficiency = 90% 

Initial cost = $ 2528 [43] 

Life time = 12 years 


Note: $1.0 = RM 3.90 


In this paper, two scenarios are investigated to determine the optimal sizing of SHRES. In the first 
case, the hybrid system consists of PV panels, wind turbines, and battery bank, while the second case 
comprised of PV panels and battery bank only. The output power curve generated by the first and second 
case of SHRES and the state of charge and discharge of the battery on an hourly basis under the best 
configuration are represented in Figure 6a and 6b, respectively. Even though the total energy generated was 
sufficient to cover the peak load in the evening, the surplus power was employed to charge the battery bank. 
These figures clearly indicated that the battery, in the period from 7 PM to 7 AM, was state of discharge and 
it satisfied the load demand, while the period from 7 AM to 7 PM was state of charge. In consonance with the 
two case studies, the contributions of output PV, wind turbine and battery bank through a one-year period are 
shown in Figure 7a and 7b. Based on Figure 6, it can be seen that the area had low wind speed and high solar 
radiation. Further, it was evident that the use of PV panels has a great advantage because the renewable 
energy for this location will enable the communities to access energy for their daily living. 



Figure 6 (a). Output of PV-WT and charge-discharge curve of the batteries 
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Figure 6 (b). Output of PV and charge-discharge curve of the batteries 



Figure 7. (a) Contribution of energy using first case study (PV, wind and battery during one year) and 
Figure. 7(b) Contribution of energy using second case study (PV and battery during one year) 


For the first case study, the set of solutions obtained from the NSGA-II and MOPSO approaches for 
one year are shown in Figure 7a and Figure 7b, respectively. Each solution represents the LPSP and COE that 
demonstrate the minimum value of the multi-objective optimization set of solutions known as a Pareto- 
optimal set or Pareto front. This technique selects one of the different solutions and makes decision on LPSP 
against COE based on the number of PV (NPV), wind turbine (NWT) and battery energy storage (NBES). 
Any of the solutions can be considered optimum, which means that no improvement can be achieved on one 
of the objective functions without aggravating the other objective function. In order to make the best 
decision, a number of points on the Pareto front was selected, and then the optimal solutions were chosen 
based on the tradeoff between cost and reliability. 

Figure 8c shows the comparison of space of operating points (non-dominated); the matching 
between two methods of optimization results was quite close. The NSGA-II algorithm had good global 
search ability but a slower convergence speed than MOPSO. 

Figure 9a and Figure 9b show the Pareto optimum final value for the second case study. Once again, 
it proved that, compared to MOPSO, the NSGA-II algorithm had a robust search capability but a slower 
convergence speed. 
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Figure. 8 (a) Pareto front for first case study using NSGA-II 



Figure. 8 (b) Pareto front for first case study using MOPSO 



LPSP(%) 


Figure 8 (c). The matching between the NSGA-II and MOPSO for the first case study 
according to pareto optimum 



Figure 9 (a). Pareto front second case study using NSGA-II 
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Figure 9 (b). Pareto front second scenario using MOPSO 


0.12 -| 

o.i 

0.08 
J 0.06 

uT 

8 0.04 


A NSGA-II 
X MOPSO 


\ 


*XA 


* A 




Optimum value of 
both methods 


1 1.5 

LPSP(%) 


Figure 9 (c). The matching between the NSGA-II and MOPSO for the second case study 

according to pareto optimum 


The optimal sizing solution obtained by the multi-objective NSGA-II and MOPSO algorithms for 
the first and second case study is illustrated in Table 3; the two case studies were compared according to the 
two algorithms, in terms of LPSP, COE and the total cost (Ctotal). When the two hybrid system designs with 
minimum LPSP and COE were considered, it was clear that the contribution from PV energy was high. 
However, the energy mix in the first case study was quite low for wind turbine compared to PV. 


Table 3. Optimal sizing result for two case studies 


Case Study 

Method 

NPV 

NWT 

NBES 

LPSP(%) 

COE ($/kWh) 

Total cost ($) 

First 

NSGA II 

178 

1 

1 

0.0856 

0.1099 

66590 


MOPSO 

178 

1 

1 

0.084 

0.1083 

65212 

Second 

NSGA-II 

177 

- 

1 

0.0998 

0.0974 

60042 


MOPSO 

178 

- 

1 

0.0966 

0.0971 

58993 


5. CONCLUSION 

This paper presented the optimal sizing of SHRES for a rural area settlement in Malaysia under the 
prevalent solar and wind conditions. The NSGA-II and MOPSO optimization techniques were used to 
optimize the sizing. The hybrid power system modeling was conducted in MATLAB using two case studies 
to select the best configuration in order to find the optimum number of variables; the first case study 
consisted of PV/WT and BES, while the second case study consisted of PV/BES to minimize 
LPSP and COE. 

The results showed that the comparison between the two algorithms was successful in supplying 
household load with two operation scenarios, and in varying the availability of solar and wind power during 
the entire one-year period. Also, the proposed performance model for the two modes were presented, where 
the NSGA-II method has more robust performance, while the MOPSO was faster in implementing the 
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optimization. The difference between these two methods in terms of number of variables, as well as LPSP 
and COE is not significant. Thus both these methods can be used for techno-economic optimization of 
SPIRES. The second case study has proved to be a more efficient solution to meet the energy demand of the 
remote areas in this location. For the wind turbines have proved to be ineffective and it's not reliable due to 
low wind speed and also high installation and operational cost, but its utilization may be highly expanded in 
the future. As conclusion, the proposed system can fulfill the targeted constraints for rural electrification and 
ensure that the load demand is met. 
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Nomenclature 


AD 

autonomy days 

Cl 

cognitive parameter 

C2 

social parameter 

Cb 

battery capacity 

DOD 

depth of discharge and equal 

Eb 

state of battery (charge and discharge) 

Ebmax 

battery maximum charge quantity 

Ebmin 

battery minimum charge quantity 

El 

hourly load demand (KWh) 

G 

global solar radiation (Kw/m 2 ) 

ir 

number of times replacement of battery, inverters, and regulators 

during the lifespan of the system 

gbest 

global best position 

k t 

temperature coefficient of the PV panel and equal 3.7*10 -3 (1/C) 

N 

size of population 

Pi 

parent population 

Pbest 

global best position 

Ppv 

output power of solar PV 

PVsTC 

nominal power under reference conditions 

Pr 

rated wind power 

PwT 

output power of wind turbine 

Qi 

population of children 

rl and r2 

represent random numbers distributed uniformly between (0 and 1) 

Ri 

populations combined (parents and children) 

Ta 

surrounding temperature 

Tc 

temperature coefficient of the PV panel 

Tref 

cell temperature and equal 25 

V 

particle velocity 

v 

wind speed 

Vci 

cut in speed 

Vco 

cut out speed 

Vr 

nominal wind speed 

X 

particle position in iteration t 

gbat 

efficiency of battery 

ginv 

inverter efficiency 

AEnet 

net energy of SHRES 

G 

self-discharge rate 
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